Boosting the Rotational Sensitivity of Matter-wave Interferometry with Nonlinearity 
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We propose a mechanism to use nonlinearity arising from inter-particle interactions to significantly 
enhance rotation sensitivity of matter-wave interferometers. The method relies on modifying Sagnac 
interferometers by introducing a weak circular lattice potential that couples modes with opposite 
orbital angular momenta (OAM). The primary observable comprises of the modal population dis- 
tributions measured at particular times. This provides an alternate mechanism for rotation sensing 
that requires substantially smaller ring size, even in the linear non-interacting regime. Nonlinearity 
can improve the sensitivity, as well as operation timescales, by several orders of magnitude. 
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Introduction: Atom interferometry has been pushing 
the hmits of precision rotation sensing in recent years [l|- 
|3|. A common perception in the field is that the inher- 
ent atom-atom interaction, of the degenerate gases used 
as the working medium, is undesirable as being another 
source of systematic error. In this paper we show that 
the nonhnearity arising from such interactions can ac- 
tually be utilized to significantly enhance the sensitivity 
of matter-wave rotation sensors. This work is motivated 
by a recent work on matter- wave transport [^, where 
nonlinearity was found to have a similar effect in linear 
geometry. 

The primary operational principle in rotation sensors, 
atomic or optical, is the Sagnac effect. The passive ver- 
sion [5| favored in atomic implementations measures the 
phase shift arising from the differential travel-time of 
counter-propagating waves in a rotating loop; while the 
active version, used in sensitive ring laser gyros [y], uti- 
lizes the frequency splitting induced by rotation in op- 
positely travelling modes in a ring-shaped laser cavity. 
Sagnac-based sensors, despite their long history, have 
seen scant considerations of positive utility of nonlin- 
earity, with a few notable exceptions for optical real- 
izations, such as a passive scheme proposal [7| to use 
non-reciprocity of nonlinear susceptibilities for interfer- 
ing beams of different amplitudes; and more recently 
there have been some advances [8|,|9[ in increasing Sagnac 
sensitivity by manipulating the anomalous dispersion of 
light. But such considerations have not found extensions 
to matter waves. 

Although it has been generally challenging to use the 
Kerr-type nonlinearity [6|, typical of both optical and 
cold atomic sensors, to enhance rotation sensitivity in 
the standard Sagnac configuration, we find a way to do 
so by introducing a crucial modification and choosing an 
alternate observable, distinct from the phase shift and 
the frequency splitting. There is, however, an indirect re- 
liance on the frequency-splitting, in which it differs from 
current matter-wave rotation sensors that rely on phase 



FIG. 1: (Color online) The atoms are trapped in a toroidal ge- 
ometry with effectively one-dimensional dynamics about the 
major circle, (a) During operation a circular lattice poten- 
tial exists as shown, (b) During readout, the lattice is turned 
off, and the counter-propagating modes are out-coupled to 
opposite directions of a tangentially coupled waveguide. 



shifts. Both attractive and repulsive non-linearities can 
be utilized. 

Model and analysis: We consider a Bose-Einstein con- 
densates (BEC) in a toroidal trap as shown in Fig. [1] 
Such a potential can be realized for example with orthog- 
onal intersection of a Laguerre-Gaussian laser beam [10[ 
with the waist of a Gaussian beam, both red-detuned 
with respect to the relevant atomic transition. Such 
traps can sustain long-lived '^ 40 s condensate super- 
fluid flow due to excitations of angular momentum eigen- 
states Uj. Assuming the minor radius of the torus to 
be much smaller than the major radius R, as is possible 
with a higher order Laguerre-Gaussian beam 1^ |l3| and 



tightly focussed Gaussian beam, the dynamical descrip- 
tion is taken to be one-dimensional (ID) along the major 
circle with excitations transverse to it frozen. The time 
evolution of the system with N condensed atoms is then 
described by the ID Gross-Pitaevskii equation in the in- 
terval X € [0, 27ri?), with periodic boundary conditions 
*(a; = 0,t) = *(a; = 27ri?,t). 



ihdt'i'{x,t) = 



2m 



+ U{x,t)+g\^{x,t)f 



^ix,t), (1) 



where m is the atomic mass and the nonlinear constant 
g ~ NgiD is defined by the number of atoms A^, and the 



effective ID atom-atom interaction strength guj [14 



For rotational sensing, the potential is taken to be of 
the form U{x,t) = U{x ~ Rilt), which describes a sta- 
tionary U{x) rotating around the center of the circle with 
angular frequency fl. Transforming to polar coordinates 
in a rotating reference frame by setting a; — >■ j/ = x — Rflt, 
y e [0, 27ri?), Eq. ([T]) becomes 



ihdt^{y,t) 



(2) 



— ^ + ihRndy 

2m 



Uiy)+g\^iy,t)\' 



*(y,0, 



with the new term ihRfldy = —flL, L being the angular 
momentum operator about the symmetry axis of the ring. 

The Raman techniques of transferring quantized or- 
bital angular momentum (0AM) from light beams to 
condensates |15| . allow one to create circulating con- 
densate with the wavefunction of the form {y\±l) = 
l/V2^exp(±i/y/i?), I e N. When U{y) = and O = 0, 
states \ztl) are degenerate 0AM eigenstates of Eq. ([2]). 
This leads to the choice of the potential to be 



U{y) ^Uo cos' (ly/R), 



(3) 



a circular optical lattice [16| (the circular analog of op- 
tical Bragg gratings widely used to manipulate linear 
atomic momentum [l7[). Such a potential can be cre- 
ated by interfering two counter-propagating Laguerre- 
Gaussian beams with opposite orbital angular momen- 
tum ±hl per photon, possible in our case by partial retro- 
reflection of the trapping Laguerre-Gaussian beam. 

Two mode approximation: We work in the limit of 
weak potential Uo <C h'P/2mR', where the atom dy- 
namics is effectively confined to a two dimensional sub- 
space spanned by the states |±/), strongly coupled with 
each other, but negligibly coupled to other states. Such a 
choice of optical potential results in the simplified atomic 
dynamics of an effective two-level system described by 



\^it))^a{t)\l)+b{t)H). 



(4) 



Inserting this into Eq. ([2]) and disregarding rapidly oscil- 
lating terms and global phase factors leads to a pair of 
coupled non-linear equations for the amplitudes: 



id(r) 

ib{T) 



ua{T) 



1_ 
ljI 



-Lol-^\a{T)\'\a{T)+uh{T), 



J_ 

2tt 



Hr)\' b{r), 



(5) 



written in dimensionless form with redefined quantities: 



T = th/{mR'), 
u = Mo/4 = UnmR^/Ah'^, 



7 EE gmR/r?. (6) 



The 'dot' signifies derivative with respect to t. Our nu- 
merical results are obtained by solving the full Eq. ^ 
which confirm that the two mode-model is quite accu- 
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FIG. 2: (Color online) Time evolution of the populations 
|a(r)p, (|6(t)P) in states \±l) and their sum, is plotted by 
numerically solving Eq. ([2]). The constancy of the sum justi- 
fies the two-mode model Eqs. ([5]) t/o <C h^l^/2mR^ , assumed 
here; breakdown occurs only at much higher lattice potentials 
(2000 times higher in (b)). With nonlinearity, there can be 
both (c) oscillatory and (d) self-trapped regimes. In all plots, 
/ = 5 and i? = lO'^m. 



rate, as indicated by the small deviation from the nor- 
malization |a(T)p -I- |6(t)P ~ 1 seen in Fig. [5] for weak 

Uiy). 

Linear regime, j — 0: When the atoms are non- 
interacting, which is possible by using Feshbach reso- 
nance to tune 7 — >■ 0, Eqs. ([5]) reduce to a pair of lin- 
ear equations which can be exactly solved. Defining 
?; = \/u'^ + uj'^P, the solution can be written as 

a{T) — [cos(77t) + i — sin(7]T)]a(0) — i— sm{riT)b{0) (7) 

6(r) — — i— sin(?7T)a(0) -I- [cos(?/r) — i — sm{riT)]b{0) . 

As seen from Eqs. ([5|), the lattice potential U{y) leads 
to transitions between the states \l) and \—l), while the 
angular velocity u lifts their degeneracy. We choose our 
initial state to be a(0) = 1,&(0) = 0. If there is no 
rotation, uj ~ 0, Eq. ([7]) implies a complete population 
swap between the states at time r^ ~ 7r/2u, i.e. a{Ts,Lu ~ 
0) = 0, b{Ts,u} = 0) = 1, as seen in Fig. [2Ia). But, for 
a; 7^ the initial state remains partially populated: 
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due to the lifted degeneracy of states |±/) in the rotating 
frame. The dependence of |a(Ts)p on the rate of rotation, 
w, is plotted in Fig. [3lja) for the same parameters used 
in Fig. [2ja). The clear dependence on uj, indicates that 
the occupation probability after a fixed evolution time, 
|a(r5)p can used for rotation sensing. 

Nonlinear regime, 7^0: The nonlinearity leads 




FIG. 3: (Color online) With |a(0)|^ = 1, the population 
|a(Ts)p of state |/) at a fixed evolution time Ts depends on 
the angular velocity u; as shown, (a) Nonlinearity sharpens 
the dependence on the u about a critical value, about u — 0. 
(b) Details of the nonlinear profiles from (a) show that the 
sensitivity to uj is over two orders of magnitude higher than 
for the linear case (note the Lo-axes scales). Time scales of op- 
eration are also shorter by a similar factor. In all plots, I — 5 
and R = 10~^m; parameters shown in (b) apply to nonlinear 
plots in (a). 



to an additional timc-dcpcndcnt potential due to the 
atomic density pattern which, in the two-mode picture, 
is ex [ab* cxp(i2ly / R) + c.c.]. Since the nonlinear terms 
in Eqs. ([5]) appear diagonally along with the rotational 
term, they directly influence the sensitivity to the rota- 
tion. General insight into the behavior can be deduced 
by noting that these equations map to those for a non- 
rigid nonlinear pendulum [18|; thus, the dynamics dis- 
plays both oscillatory and self trapped regimes as shown 
in Fig. [2Uc) and (d) . Nonlinear strength above a critical 
point (7 > 2ttuq for 7 > 0, w = 0) causes self trapping 
that limits the transitions between states |±^) resulting 
in only partial transfer during a cycle; eventually there 
is no transfer at all for sufficiently high nonlinearity. 

The nonlinearity can be used to significantly enhance 
rotation sensitivity by magnifying the dependence of the 
primary observable |a(Ts) p on u) as we now explain. For a 
fixed potential u, and for w = 0, the nonlinearity is tuned 
to be just below the critical point in the oscillatory regime 
and the time required for complete swap Ts is set to be the 
detection time. Now, if there is rotation uj ^ 0, there are 
two possibilities depending on the direction of rotation: 
(i) In one direction, the rotation tips system over into the 
self trapped-regime so that the |a(rs)p increases from 



for complete swap at w = to 1 for no swap at all for a 
tiny change in uj; (ii) When the rotation is in the opposite 
direction, the system stays in the oscillatory regime, but 
due to the proximity to the critical point, there is still a 
high sensitivity to small changes in uj. This behavior can 
be seen in Fig. [3l Changing the sign of the nonlinearity 
switches the dependence on the direction of rotation in 
the relevant regime close to cj = 0, which follows from 
the relative signs of w and 7 in Eqs. ([5|). 

The most remarkable feature apparent in the figures is 
that compared to the linear case, the sensitivity (for the 
assumed nonlinearity) is enhanced by over two orders of 
magnitude. An additional attractive feature implicit in 



the parameters shown in the figures, is that for experi- 
mentally feasible ll| values of 7, the associated values of 
the lattice potential u of interest are substantially higher 
than in the linear case; this has the effect of reducing the 
swap time Tg by about two orders of magnitude as well. 
Thus, nonlinearity not only increases the sensitivity, but 
the experiment can be accomplished in shorter times as 
well. 

Rotational Sensitivity: The complementary probabili- 
ties of the two modes indicates a binomial (or Poisson) 
distribution, with each of the TV atoms constituting an in- 
dependent trial. This can be brought about by tuning the 
nonlinearity to zero using Feshbach resonance during the 
detection process. Therefore, a criterion for sensitivity 
to rotational frequency w is provided by comparing the 
binomial mean value N\a(Ts)\'^ and the standard devia- 
tion v]V|a(Ts)&(Ts)| for u! = 0, with that of the minimum 
value u! = Aw that is distinguishable from zero. This 
sets the criterion for frequency resolution and sensitivity 



A^laol 



N\aobo\ < N\aAu 



N\aAujbAu. 



(9) 



where a^, hi are the respective amplitudes of the modes. 
Using this criterion with parameters from the experi- 
ment in Ref. ll| with sodium atoms, m = 3.8 x 10~^^ 



kg, R 



10- 



N 



2 X 10 , for our choice of 



/ = 5, we compute the short term sensitivity in the 
linear case to be Aw = 1.1 x 10^"'' which corresponds 
to AQ = 2.9 x 10~* rad/s, with a run-time tg = 7.6 
s. In the nonlinear case, the experimental parameters 
lead to 7 = 96 as used here, and the sensitivity is sig- 
nificantly enhanced Aw = 2.5 x 10~^ which corresponds 
to Afl = 6.5 X 10^^ rad/s, and the run-time shortened 
commensurately tg — 0.087s. By repeating the cycles 
we can define a process with cumulative atom-number 
and duration, so that on normalizing by the atoms used 
per second, Afl x \/t^, we can frame our sensitivities 
in commonly used units: For the linear case 2.8 x lO^"' 
rad/s/-\/lIz and for nonlinear 1.9 x 10^^ rad/s/\/lIz. 

In principle, by choosing the nonlinearity to be increas- 
ingly closer to the critical point, the sensitivity can be in- 
creased arbitrarily by causing more pronounced changes 
to |a(Ts)p about w = 0. The limiting factor is the con- 
current heightened sensitivity to the potential u. This 
means that any fluctuations in u will change the position 
of the dip in Fig. [3l so that for the same r^ the mini- 
mum of the curves would be shifted Suj from w = 0. In 
Fig. |4]we present a log-log plot of 6u/u versus Suj, deter- 
mined by numerical time-evolution for u ± Su, and find 
that the dependence is essentially linear. In cold atom 
experiments stability of {Su)/u ~ 10^^ can be reached. 
To access the enhanced sensitivity attainable with non- 
linearity, much better stability would be needed, such as 
attainable in some highly sensitive gravitational experi- 
ments, where stability of {5u)/u ^ 10~^ is within reach 

M. 




FIG. 4: (Color online) Phase space diagram of probability dif- 
ference versus phase difference of the mode amplitudes a{Ts) 
and 6(ts). Each line represents time evolution for specific val- 
ues of 7 and uj. (a) The inset shows transition from oscillatory 
to self-trapped behavior as interaction strength varies through 
the critical point at 7 ~ 96.1 for fixed uj = 0. With the inter- 
action strength fixed at 7 = 96 instead (just below the critical 
point for lu = 0) in the primary plot, similar behavior occurs 
when w is varied, (b) Deviation Sco, of the location of the 
minimum in Fig. [3l from a; = is linearly proportional to 
the stability 5u/u of the lattice strength, the linear fit being 
log(5aj) = 1.04 X log{Su/u) + 0.011. 



Comparison to the Sagnac effect: If both the nonlin- 
earity and lattice are absent, then Eq. ([5]) shows that 
rotation simply splits the energies of | ± Z) levels by 2wZ, 
and we find that the corresponding frequency splitting in 
physical units is 5v = [AViA) / {\P) which is identically 
the Sagnac frequency splitting for matter waves, with 
A = TrR^,P = 2ttR the area and perimeter of the ring 
and A = h/{mv) the de Broglic wavelength. This implies 
that the physical process in this limit is identical to the 
Sagnac effect. Even without rotation a; = 0, the nonlin- 
earity causes an energy splitting (7/27r)(|a(0)p — |6(0)p), 
(the populations remain constant due to lack of coupling 
in the absence of the lattice). This causes the standing 
wave formed by the counter-propagating modes to rotate. 
Any physical rotation w 7^ also causes a similar rota- 
tion which will simply be added on to the nonlinearity- 
induced rotation. Since the nonlinearity itself does not 
depend on the rotation, it does not help in rotation sens- 
ing and is in fact undesirable because now the effect of the 
actual rotation has to be distinguished form a potentially 
much stronger fringe rotation due to the nonlinearity. 

This makes it clear that the coupling introduced by the 
lattice is the crucial element in our proposed mechanism. 
It leads to the population oscillation that provides the 
alternative observable in the population of the modes, 
distinct from the phase shift or the frequency splitting 
that has been traditionally used in Sagnac-based rota- 
tion sensing. Crucially the coupling makes it possible for 
nonlinearity to have an enhancing effect on the sensitiv- 
ity, which has not been possible in the traditional Sagnac 
setup. 

Detection of the motional states: Rotational sensing in 



our scheme requires precise determination of the popu- 
lations of states |±/). One possible way is indicated in 
Fig. [Ijb) that can work as follows: After duration r^ of 
switching on, the lattice potential will be turned off and 
the nonlinearity 7 tuned to zero. This freezes the popu- 
lations in the respective rotational states and avoids any 
complications due to nonlinearity during measurement. 
Then, a red-detuned Gaussian beam is coupled tangen- 
tially to the toroidal trap allowing the atoms to tunnel. 
The atoms in the two states | ± I) will move out in oppo- 
site direction along the tangential beam and can be sep- 
arately detected to yield the ratio of |a(r)p and |6(r)p. 

Conclusions: We have proposed a way to use non- 
linearity to boost rotational sensitivity of matter-wave 
interferometers in a configuration where the observable 
comprises of populations of coupled modes at the end 
of a operation cycle of fixed time. This provides a way 
around the challenges faced in attempts to use nonlin- 
earity to improve sensitivity of Sagnac interferometers. 
Even in the linear regime, our method offers high sen- 
sitivity levels that could be substantially improved with 
increasing particle number. Nonlinearity can boost the 
sensitivity to levels competitive with current matter- wave 
interferometers |20[, and notably this is achieved with a 
ring of area ^ 10~^ mm^ that is smaller by several orders 
of magnitude [2, l3| . The nonlinear enhancement requires 
heightened stability of the mode-coupling potential, and 
can have applications impact in highly sensitive experi- 
ments where such stability is accessible. 
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